Charge-dependent migration pathways for the Ga vacancy in GaAs 
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Using SIEST-ART, a combination of the local-basis ab-initio program SIESTA and the activation- 
relaxation technique (ART nouveau) we study the diffusion mechanisms of the gallium vacancy in 
GaAs. Vacancies are found to diffuse to the second neighbor using two different mechanisms, as well 
as to the first and fourth neighbors following various mechanisms. We find that the height of the 
energy barrier is sensitive to the Fermi-level and generally increases with the charge state. Migration 
pathways themselves can be strongly charge-dependent and may appear or disappear as a function 
of the charge state. These differences in transition state and migration barrier are explained by the 
charge transfer that takes place during the vacancy migration. 

PACS numbers: 61.72.Ji, 71. 15. Mb, 71.55.Eq, 71.20.Nr 
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I. INTRODUCTION 

Self-difFusion is one of the basic mass-transport mech- 
anisms in materials. While it is one of the most power- 
ful tools used in the preparation of nanostructuresiiSiSii, 
many questions remain regarding the microscopic de- 
tails of self-diffusion and recent studies have shown that 
even in the simplest cases, the mechanisms can be much 
more complicated than was initially thought. Diffusion 
on simple metallic surfaces, for example, was found to 
take place through a range of mechanisms involving from 
one to at least seven atoma^. Similarly, recent studies 
on self-interstitial clustersL^ and disordered systems2ii£ 
have highlighted the importance of collective moves in 
easing atomic motion even in bulk systems and under- 
lined the importance for a detailed characterization of 
these mechanisms in materials of technological interests. 
This is the case for semiconductors, for example, that are 
at the heart of the electronic industry. While one pre- 
dicts, using symmetric considerations, the self-diffusion 
pathways in elemental materials, such as silicon, this ap- 
proach becomes impossible when impurities are added or 
multi-compound systems are considered. Thus, one must 
turn to experiments or numerical simulations to provide a 
direct and comprehensive study of diffusion mechanisms 
in semiconductors. The difficulty to extract precise in- 
formation on the diffusion mechanisms in these materials 
is compounded by the nature of bonding and the impor- 
tance of charged defects that complicate seriously both 
experimental and theoretical studies. 

Previous numerical studies of the migration pathways 
of intrinsic defects in binary semiconductors have focused 
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on GaAsii, SiC— ii^ and GaNi4. Except in one caseii, 

where high-temperature molecular dynamics was used, 
these works focus on optimizing preselected pathways us- 
ing algorithms such as the dragii or the ridge methodic 
that rely on the knowledge of the initial and final states in 
addition to a decent guess of the overall diffusion mecha- 
nism. While these approaches work efficiently to identify 
with precision the migration energy of previously known 
diffusion trajectories, they cannot help identify complex 
or unexpected mechanisms that could also play an im- 
portant role in the diffusion process. 

Here, we present the application of the activation- 
relaxation technique in its more recent implementation 
(ART nouveau) to explore systematically the diffusion 
mechanisms of complex systems. More precisely, we fo- 
cus on the diffusion of Voa in GaAs because of its ap- 
parent simplicity but also because of its technological in- 
terest and its role in affecting the properties of bulk ma- 
terials and nanostructures. Gallium vacancies are found 
to be mobile at typical growth and annealing tempera- 
turesi^ — with a dominant charge state strongly depend- 
ing on growth conditions, temperature, dopants, etc. — 
playing the main role in dopant diffusion. For example, 
Bracht et alm^ showed recently that the contribution of 
Vca to Ga self-diffusion in GaAs is even more important 
than earlier estimation giving an important contribution 
to the total diffusion profile. Furthermore, Tuomisto 
et ali^ found that Ga vacancies play a central role in 
the migration of Mn in Gai.xMuxAs alloys. Finally, the 
energy blueshift in PL spectra of InGaAs/GaAsiS and 
InAs/GalnP quantum dotsSS has also been recently at- 
tributed to Vca diffusion. 

Focusing on a simple defect, a gallium vacancy (Vca), 
in the weakly ionic GaAs, we show here that diffusion 
in bulk semiconductors is a complex phenomenon that 
depends not only on the geometry of the defect and the 



surrounding lattice but also on its charge. In particu- 
lar, we identify a new mechanism for the diffusion to the 
second neighbor in addition to the one already found by 
Bockstedte and Scheffler— , plus a number of other jumps 
to the first and fourth neighbors. Not all these pathways 
are likely to occur in a normal range of temperature, and 
some exist only for a subset of charge states, but their 
existence underlines the under-estimated richness of dif- 
fusion mechanisms in bulk materials. 

This paper is organized as follows, in Scc^I we present 
the details of activated event generation using SIEST-A- 
RT. A description of initial and different transition states 
in various charge states are presented in Sec. IIIII In 
Sec. IIVI we discuss the origin of charge dependent mi- 
gration barrier and we compare our results with experi- 
mental and theoretical works. 



II. DETAILS OF THE SIMULATION 

Our simulations are performed using SIEST-A-RT, a 
method combining a self-consistent density functional 
method (SIESTA )2i with the most recent version of acti- 
vation relaxation technique (ART nouveau)^. Integrat- 
ing various empirical potentials, ART nouveau was shown 
to sample efficiently the energy landscape of amorphous 
semiconductors^, glassesS^, and proteins^, for example. 

SIEST-A-RT was used for the first time to study va- 
cancy diffusion in Si and details of its implementation 
can be found in Reim^. Forces and energies are eval- 
uated using density-functional theory (DFT) with the 
local-density approximation (LDA) using standard norm- 
conserving pseudo-potentials of Troullier-Martins2& fac- 
torized in the Kleiman-Baylander form^. Matrix ele- 
ments are evaluated on a 3D grid in real space. The 
one-particle problem is solved using linear combination 
of pseudo-atomic orbitals (PAO) basis set of finite range. 
Here, we use the extended local basis seliS^, which was 
shown to reproduce closely the best solution while mini- 
mizing computational costs. Calculations are performed 
on a 215-atom GaAs supercell sampled at the F spe- 
cial point. All relevant charge states (0, —1, —2 and 
—3) are fully relaxed until the residual force falls below 
0.002 eV/A , then we proceed to search for local diffu- 
sion pathways by assuming that the charge state of the 
defect is preserved during any transformation. Unless 
mentioned, all events start from fully relaxed gallium va- 
cancy geometries. 

In order to break the initial local symmetry, activated 
events are started by displacing in a random direction a 
region of the cell centered around a chosen atom in the 
first, second or fourth-shell around the vacancy, involv- 
ing typically between 10 and 30 atoms. The structure 
around the vacancy is then deformed along this random 
direction until the lowest curvature, corresponding to the 
lowest eigenvalue of the Hessian matrix, becomes nega- 
tive, falling below a preset threshold value. The system is 
then pushed along the corresponding eigenvector, while 




FIG. 1: (Color online) Constant electron density surface at 
0.0004 electron/A^ which shows the degree of localization of 
extra —3 charge around a gallium vacancy in GaAs located 
initially at the center of the box. 

minimizing the energy in the perpendicular hyperplane, 
until the total force falls below 0.1 eV/A, indicating that 
the transition state has been reached. 

About 60 events were generated in total for all charge 
states, with 20 events for jumps to the first nearest neigh- 
bor and 30 for jumps to the second neighbor. In both 
cases the structure is deformed by at least 0.9 A be- 
fore a sufficiently negative eigenvalue appears. Diffusion 
to the fourth neighbor was more difficult to complete (we 
generated 10 such events in total) and we set up a thresh- 
old displacement of 1.4 A to allow the configuration to 
escape from the harmonic well. 



III. RESULTS 
A. Gallium vacancies at the initial state 

We first examine the vacancy at the initial state before 
moving to the analysis of transition states. The struc- 
ture of Vca for all charge states, for a number of basis 
sets within SIESTA, is presented in Refi2& and agrees 
well with other simulations and experiments. Relaxed 
at its energy minimum, Vca conserves the symmetry 
for all charges, while the open volume associated with the 
vacancy decreases with increasing charge. Spin-polarized 
LDA relaxations for Vca also lead to a Td symmetry, con- 
firming that the symmetry conservation of Vq^ is not a 
drawback of LDA but a behavior proper to cation vacan- 
cies, as pointed out previously by Chadi^S. This behav- 
ior has been further confirmed by a recent calculation's 
on Voa in GaAs using two-component density-functional 
theory-LDA applied to a cubic 64 atom supercell together 
with a 4^^ k-point mesh Brillouin zone sampling . 

A detailed analysis of the charge densities plots (see 
Fig. ^ and the MuUiken populations reveals that the 
electrons added to the neutral vacancy are delocalized 
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TABLE I: Nearest neighbor distances (in A) in 111 direction relevant for first nearest neighbor diffusion from the initial — » 
to the final state . Distances are calculated by taking the initial position of Vca as reference. The last column describes the 
geometry of the final state. 

Vca - As^"* As^'' - Ga^""^ Ga^"'' - As^'''' Final geometry 



VSa 2.08^0.44 2.43^3.40 2.45 ^2.46 (Asoa + Vif)" 

V^^ 2.08^0.90 2.42^2.65 2.45^2.58 (Vca + Ias + V^f)-^ 

2.06^0.53 2.41^2.69 2.45^2.64 (^Asoa + VXf + loa + Vta''y^ 

Vq^ Does not diffuse to the first neighbor 



and distribute themselves on the 111 axes passing by the 
As dangling bonds in a similar way for all studied charge 
states. Less than 4% of the charge is localized on the four 
As neighboring the vacancy {As^^*') while the remaining 
96 % is spread over more distant neighbors on the 111 
axes. This suggests that these predominantly covalent 
As— Ga bonds arc progressively stiffened upon the addi- 
tion of extra electrons, thus making them more difficult 
to break compared to the remaining bonds and affecting 
directly the height of the diffusion barriers. 



B. Diffusion path to first neighbor 

The activated events we generated using SIEST-A-RT 
show that diffusion to the first nearest neighbor is not 
possible for all sequential charge states, contrary to what 
was proposed by Van VechtenSi nor impossible, in con- 
tradiction with what was found by Bockstedte and Schef- 
fler— . We find rather that diffusion to the first nearest 
neighbor is very much charge dependent. 

For g = the first neighbor of the vacancy (As^**) 
diffuses, along the 111 direction, toward the vacant site 
via a split vacancy configuration by optimizing its bonds 
with the close neighbors. This mechanism is similar to 
the diffusion mechanism of a neutral silicon vacancy in 
gjjj(,Qjj25j32j33^ First, each back bond of the diffusing atom 
(As^** — Ga^"'*) stretches until it breaks during the mi- 
gration of As'^"*' toward the vacancy. As^'^^ proceeds in 
its migration until a metastable vacancy-antisite struc- 
ture {AsGa + VXf) forms. The As^** - Ga^"'' bond 
evolves from 2.43 A at the initial state to 3.40 A at the 
final state. This metastable complex is 0.67 eV higher 
than Voa and the recorded migration barrier for the 
first neighbor diffusion is 0.84 eV. The next first neigh- 
bor jump [V^f Ga^"'*) leading to the formation of 
AsGa + Ga^^l + complex is also possible by cross- 

ing a barrier of 1.55 eV. Starting from this last complex 
we did not find any mechanism for Vj"'' to diffuse fur- 
ther to the first neighbor, suggesting that this jump is 
unfavorable. 

For q = —\ the vacancy follows the same path to the 
saddle point as for g = 0: it crosses a barrier of 0.9 eV 
passing by a split vacancy site, but As^** — Ga^"'' bond is 



less stretched. As^'^^ relaxes then in a split vacancy site, 
forming a (Vba + Ias + VlfY^ complex. As"'^ - Vca 
distance is reduced to half (2.08 A ^0.90 A) suggesting 
that As^"* is at half way between the two vacancies. This 
metastable configuration is 0.81 cV higher than the ini- 
tial minimum. We confirmed that this final metastable 
state is not only a local minimum along the diffusion path 
to the first neighbor by starting from an ideal vacancy- 
antisite complex {AsGa + Kil*)"^' then relaxing until 
the residual force becomes lower than 0.002 eV/A. This 
vacancy-antisite structure is found to be unstable in — 1 
charge state because As^^^ leaves the ideal antisite struc- 
ture and prefers to relax at a metastable state half way 
between the two vacancies. 

For q = —2, we find only a collective motion of a 
As^"* — Ga?'^'^ pair toward the vacancy along 111 direc- 
tion. As^^* — Ga'^"''^ bond stretches slightly, while As^"^* 
approaches Vq^ as close as 0.53 A. Consequently, Ga?'^'^ 
is forced to stretch its back bonds and stabilizes in an 
interstitial position. Finally, As^''* atom occupies Vca 
while Ga'^"'^ is located at a split interstitial position be- 
tween V^f and Vq^'^. The resulting metastable complex 
{AsGa + Vlf + iGa + Vj^'')-^ is 1.74 eV higher in en- 
ergy than Vq^ and can be obtained by crossing a barrier 
of 1.86 eV. Relevant distances for successful first nearest 
neighbor diffusion in q = 0,-1,-2 are summarized in 
Table ni 

For q = —3, all attempts for first nearest neighbor 
diffusion failed and the configuration always relaxes back 
to the initial minimum. Even when forcing the jump by 
using the transition state at the neutral charge state as 
starting point for a convergence of the Vq^ to its saddle 
point, the vacancy systematically returns to its original 
state. 



C. Diffusion path to fourth neighbor 

The situation is almost opposite for the diffusion to 
the fourth neighbour: all negatively charged vacancies 
succeed in diffusing along this pathway, only the neutral 
vacancy refuses to go this way. 

The fourth neighbor of the vacancy (Ga**'* ) approaches 
the interstitial region near the vacancy by diffusing along 
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FIG. 2: (Color online) Migration path of Vq^ to the fourth 
neighbor along (100) direction (see the text). 

the 100 direction. Figure [21iUustrates this configuration: 
the dashed arrow shows the direction of the jump from 
the initial state to the saddle point, while the full ar- 
row shows the path from the saddle to the final. Ga^*'* 
diffuses to an unstable interstitial position close to the 
vacancy. In the cubic zincblende structure of GaAs, 
with lattice constant a and a vacancy located initially 
at (0,0,0), the fourth neighbor diffuses first from (a, 0,0) 
to the interstitial position at (a/2, 0,0). The transi- 
tion state of this mechanism can be described as a gal- 
lium interstitial placed between two distant vacancies 
{Vaa - la'a ~ ^Ga)- This interstitial atom is not a di- 
rect neighbor of either of the two vacancies because the 
Qg^th _ distance shortens from 5.6 A at the initial 
minimum to 2.88 A at the transition state. This diffusion 
mechanism require an elevated barrier (around 4.24 cV) 
since the Ga^*^ needs to break two bonds initially with 
As^*^ farthest away from the vacancy and to twist the 
remaining ^.s'^'''' bonds. Once it arrives at the saddle 
point two new bonds with As^** arc formed. It is inter- 
esting to note, as is shown in Table ITTll that the energy 
barrier is almost independent of the charge for the de- 
fects that manage to diffuse to the fourth neighbour; this 
weak dependence can probably be attributed to the con- 
stant electronic charge distribution along 100 direction 
for different negative charging. 



D. Diffusion path to second neighbor 

By exploring the energy landscape of the system, we 
find that vacancy migration to the second neighbor oc- 
curs by two mechanisms. The simplest diffusion path- 
way is already well knownii and is considered to mediate 
self-diffusion in binary semiconductors, while the second 
pathway, which has not been reported, to our knowledge. 




FIG. 3: (Color online) Diffusion to second neighbor by plane- 
passing mechanism (refer to the text). 



for this system, is more complex and involves the corre- 
lated motion of many atoms neighbouring the vacancy. 



1. Plane-passing mechanism 

The most intuitive diffusion pathway, which we call the 
plane-passing mechanism, brings one Ga^"^^ to the inter- 
stitial region, joining the diffusing atom and the vacancy, 
along the 110 direction. The diffusing Ga^"'^ atom must 
go through the diffusion plane perpendicular to the 110 
direction during its way to the vacancy site causing three 
As^"^ atoms to move away from the vacancy and open- 
ing the cage. The transition state of this mechanism, 
first identified by Bockstedte and Scheffler— , can be de- 
scribed as a gallium interstitial placed between two close 
vacancies (Vca + IcTa + ^alf)- Figure 01 illustrates this 
configuration: the dashed arrow shows the direction of 
the jump from the initial state to the saddle point, while 
the full arrow shows the path from the saddle to the final; 
the diffusion plane formed by the second nearest neigh- 
bors of the vacancy is also shown. 

From Tabled it is clear that, as with the other mecha- 
nisms to first and fourth neighbors, the structural details 
of the jumps are charge-dependent. The position of the 
transition state, in particular, changes by 20 % as the 
vacancy goes from neutral to a charge of —3. As can 
be expected, the displacement of the unstable interstitial 
position as the diffusing atom becomes a nearest neigh- 
bour of the vacancy has noticeable impact on the energy 
barrier, which goes from 1.7 cV for neutral and —1 charge 
state to 1.84 and 2.0 eV for —2 and —3 charges respec- 
tively. 

More precisely, for Vq^ and Vq^ ^ the saddle point is 
a site close the hexagonal interstitial, while it is closer 
to the tctrahcdral interstitial configuration for the more 
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FIG. 4: (a) Migration trajectory to the second neighbor by 
the plane-passing mechanism for —3 charge state. 11 images 
(circles) are relaxed with CI-NEB method until a force toler- 
ance of 1 eV/A is reached, lines are guide to the eye. The 
contributions from different moving atoms are decomposed in 
two: (b) The path followed by the diffusing atom; (c) Ar- 
rangement of the other atoms around the defect. 



negative Vq^ and Vq^. At the transition state, as the 
charge increases, the displacement of the gallium atom 
out of the plane becomes less pronounced; for the —3 
charge state, it almost vanishes, leaving the moving atom 
on the plane. 

This can be seen by looking at the the full migration 
trajectory for —3 charge state in Figure 21 This path 
is generated by initially interpolating between the ART- 
generated initial, saddle and final states, generating 11 
images. These image configurations are then relaxed us- 
ing the Cl-nudgcd-elastic-band method^ until the total 
force becomes lower than leV/A . After decomposing the 
total path into contributions coming from different mov- 
ing atoms one can notice that the path followed by the 
diffusing atom is nearly symmetric, while most asymme- 
try in the total path comes from arrangement of the other 
atoms around the defect. 

Differences in the migration barriers arc mainly due to 
the diffusion of the electronic charge during the jump. 
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FIG. 5: (Color online) Mulliken population analysis of charge 
distribution around the vacancy at initial (square) and the 
saddle point (circle) of plane-passing mechanism. Shown the 
degree of localization of extra electrons at typical distances 
from the vacancy for charge —1 (top panel), and —3 (low 
panel). The insets show the constant electron density sur- 
face at 0.00018 and 0.0008 electron/ near V^^ and V^';^ 
respectively. 



TABLE II: Evolution of the distance (in A) between the initial 
position of the vacancy and the diffusing atom (Vq^ — Ga^^"^) 
in the initial state and the transition state for both diffusion 
mechanisms to the second neighbour. 





Plane- 


passing 


Cluster-assisted 




Initial 


Saddle 


Initial 


Saddle 


vL 


3.89 


1.89 






Vaa 


3.88 


2.09 


3.88 


2.5 


Vaa 


3.82 


2.19 


3.82 


2.78 




3.80 


2.35 


3.80 


2.89 



For the singly negative vacancy, Mulliken population 
analysis at the saddle point (top panel in Figure [SJ is 
compared to the initial configuration. It shows the dif- 
fusion of the electronic charge during the migration of 
the Ga^"'' atom. The electronic charge at Ga^"'' diffuses 
with it to the hexagonal interstitial site thus saturating 
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partially the As-'^*'* dangling bonds. The charge is conse- 
quently suppressed from the neighborhood of the diffus- 
ing atom (Iq^) and spreads uniformly over more distant 
shells as can be seen from the 3D charge density (top 
inset of Figure 13 ■ 

For Vq^'~^ the diffusing atom is less engaged toward 
the vacancy: circles in the lower panel in Figure |S1 are 
shifted to the right as a signature of volume opening 
that affects the second neighbor of the vacancy as well. 
However, the 3D charge densities for Vq^ and V^^ are 
different from Vq^. After the jump, the charge becomes 
highly localized around the dangling bonds belonging to 
As atoms farthest away from Iq^'^ labeled Asyj^^. Some 

of the charge carried by two closer Aslf^^^ to get 
transferred to the other Asy*^; these have non saturated 
bonds that arc still pending. 

2. Cluster-assisted mechanism 

Negatively charged Vca can also diffuse on the Ga sub- 
lattice through a mechanism that we call the cluster- 
assisted path. As far as we know, this diffusion pathway 
had not been reported until now. Instead of crossing 
directly the diffusion plane, Ga^"'^ approaches the inter- 
stitial region far from the plane, being assisted by two 
As atoms — respectively first (As^^*) and third neighbors 
{As^^'^) of the vacancy — and one gallium atom second 
neighbor of the vacancy Ga^^"^. From the initial to the 
transition state (dashed arrows on Figure |SJ), an incom- 
plete bond exchange mechanism of type Wooten-Winer- 
Weaire^^ — which we also find in Si — occurs between 
Ga^"'' and As^^'', then Ga^''^'^ is pushed into the intersti- 
tial region. 

The cluster formed by As^'* + I^^ + As^'''^ + Ga^"'^ 
plays the major role for diffusion since the bond distances 
& - ^s^""* and I^l"^ - As^'''* remain unchanged ( 2.4 A) 
regardless of the charge state, while the Vq^ — Iq^^ dis- 
tance increases by adding extra electrons (see Table . 
This assumes that the whole cluster becomes less engaged 
toward the vacancy when passing from —1 to —3 charge 
states. During the relaxation from the saddle point to the 
final state (full arrow on Figure 0) only the gallium atom 
at the interstitial position {I^^) continues its motion to- 
ward the vacancy leaving the remaining constituents of 
the cluster close to their initial positions. This is con- 
firmed by looking at the full migration trajectory for 
—3 charge state for the cluster-assisted mechanism plot- 
ted in Figure [T] From the initial to the saddle point 
many atoms are experiencing rearrangements and dis- 
placements, but the main contribution comes from the 
diffusing atom which is less engaged toward the vacancy 
than for plane-passing mechanism, while during the re- 
laxation from the saddle point to the final state the main 
contribution to the total path comes from the relaxation 
of I'^^ atom (Figure Hb)). 

MuUiken population analysis shows that the electronic 
charge diffusion for cluster-assisted mechanism is com- 




FIG. 6: (Color online) Diffusion to second neighbor by 
cluster-assisted mechanism (refer to the text). 



TABLE III: Calculated diffusion barriers (in eV) for 
Vq^^'^'^''^'^ in GaAs for all possible migration paths iden- 
tified. The empty cell means that the migration was not pos- 
sible with this mechanism. 





1/0 


Vaa 


l/~2 




First neighbor 


0.84 


0.90 


1.86 




Plane-passing 


1.7 


1.7 


1.85 


2.0 


Cluster-assisted 




2.44 


2.89 


3.24 


Fourth neighbor 




4.24 


4.24 


4.3 



pletely different from the plane-passing mechanism. By 
increasing the charge state of the vacancy from —1 to 
—3, positions of the cluster atoms are all shifted to the 
right confirming that the whole cluster is experiencing a 
displacement far away from the vacancy. In addition, the 
electronic charge becomes more and more concentrated 
around the cluster as the charge on the defect increases. 
Consequently, the migration barrier increases in signif- 
icant way with the charge state reaching 2.44, 2.89 to 
3.49 cV for —1, —2, —3 charges respectively. 

No cluster-assisted path was found at the neutral state: 
we tried to generate this event by starting directly from 
the transition state at —1 charge, but, the neutral va- 
cancy relaxes back to the initial minimum rather than 
diffuse to the second neighbor site, demonstrating that 
this path is impossible at the neutral charge state. 



E. Summary of the results 

The summary of the calculated migration barriers is 
presented in Table IIIII First neighbor diffusion has the 
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FIG. 7: (a) Migration trajectory to the second neighbor by 
the cluster-assisted mechanism for —3 charge state. 11 im- 
ages (circles) are relaxed with CI-NEB method until a force 
tolerance of 1 eV/A is reached, lines are guide to the eye. 
The contributions from different moving atoms are decom- 
posed in two: (b) The path followed by the diffusing atom; 
(c) Arrangement of the other atoms around the defect. 



lowest barrier for and —1 charges, but cannot be consid- 
ered to be dominant for Vca since the complete diffusion 
process is impossible to achieve even for the neutral va- 
cancy. However, this mechanism induces the formation 
of defects belonging to the Asoa family of great impor- 
tance since they are responsible of EL type defects in 
GaAs (see the discussion below). 

Self-diffusion will rather be dominated by jumps to the 
same sublattice. Gallium vacancies can diffuse to the sec- 
ond neighbour through the plane-passing mechanism for 
all charge states by crossing barriers lower than 2 cV. The 
impact of the charge state on the trajectory is reflected 
on the charge-dependent barrier-height energy, for exam- 
ple, the barrier height for Vq^ is 18 % higher than for 
Vq^, making the crossing 30 times less likely at 1000 K. 
The simplicity of this path as well as the relatively low 
migration barrier make the plane-passing mechanism the 
potential candidate to mediate vacancy self-diffusion in 
GaAs. 



Negatively charged vacancies might also follow the 
cluster-assisted pathway to diffuse to the second neigh- 
bor or diffuse directly to the fourth neighbor with energy- 
barrier which is, by about 50 % and 100 % respectively, 
higher than that of the plane-passing mechanism. Al- 
though these two mechanisms have elevated barriers, 
their existence is very interesting by itself as it shows 
the underestimated richness of self-diffusion phenomena 
in bulk semiconductors. 



IV. DISCUSSION AND COMPARISONS 
A. Self-diffusion to first-neighboring sites 

In a binary lattice, vacancy hops to the first-neighbor 
sites require a complex sequence of moves to preserve 
chemical ordering on the long run. In 1984, Van 
VechtenSi proposed a model leading to vacancy diffusion 
to the second neighbor in GaAs. His model describes 
a 11 first-neighbor hop process on a six membered ring 
and is based on two assumptions: (i) jumps to the first 
neighbors are always possible, and (ii) vacancy-antisite 
complexes are always stable. During successive hops the 
vacancy should proceed leaving behind a chain of anti- 
sites, which is energetically unfavorable. This is avoided 
by the vacancy passing twice through the same sixfold 
ring. At the final stage all the antisite defects are re- 
moved. Otherwise such a mechanism would create an 
unfavorable excess of antisite defects beyond the equilib- 
rium concentration. Bockstcdte and Schcfflcrii- studied 
the validity of these assumptions — possible diffusion via 
first neighbor hops — using the drag method and LDA. 
They found that Asoa + V\f is metastable in the neu- 
tral state while it is unstable for other negative charge 
states for a 64-GaAs atoms system. They concluded that 
diffusion by first nearest neighbor hops was impossible. 

Our study shows that the complex is rather metastable 
in —1 charge state too. The first nearest neighbor diffu- 
sion of Vca in —1 and —2 charges is found to be possible 
via deformed structures belonging to Asca family. While 
Vca diffusion to the first neighbor lead to a symmetric 
AsQa ~ Vas complex in the neutral state, this configura- 
tion is distorted in the case of —1 charge states. In the 
deformed structure for —1 charge state, As atom do not 
occupy the tetrahedral gallium vacancy site, it is rather 
is displaced from the Ga vacancy site by about 0.9 eV 
toward V\f. Interestingly, this deformed structure is no 
longer metastable for —2, the transition state originates 
from a AsQa + GaAs + Vca complex, whith GaAs pushed 
toward Vca and occupying a split vacancy site. 

We attribute this charge-dependent first-neighbor dif- 
fusion to two competing factors: (1) the progressive in- 
crease in the strength of As^"*" — Ga'^^'^ bonds as electrons 
are added make them more difficult to break, and (2) the 
electron density in the region surrounding the vacancy — 
especially on As dangling bonds — increases as electrons 
are added to the relaxed system. Thus, Coulomb repul- 
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sion between three As""^*** atoms and the diffusing As^'^* 
atom can become so strong that this atom cannot ap- 
proach further the vacant site. This picture causes the 
As^'^* atom to relax at the spHt vacancy configuration 
in —1 charge state. If the system is charged —2, the 
Afii"* - Ga^"'' bond is so strong that As^"* puUs Ga^"'' 
with it during its diffusion to the saddle point position 
leaving behind V^f and Vq^'^. The electronic density 
is partially transferred from the initial vacancy to these 
two new vacancies allowing As^^* to relax on the initial 
vacant site and leaving the Ga^"'^ atom stacked between 
Vif and . 

Similarly deformed structures for negatively charged 
defects have been recently reported in the litterature. 
Chadi^ found that the isolated As antisite structure 
{AsQa is generally accepted to be the basic structure of 
EL2 defect) exists in charge states +2, +1 and occupy- 
ing the tetrahedrally symmetric position and in —1,-2 
charge states when it undergoes a small displacement 
that causes a deviation from symmetry. The main 
difference between different charge states is the degree of 
relaxation that the antisite and its direct neighbor un- 
dergo. This possible high negative charging of Asca hi- 
ducing a structural relaxation is similar to the deformed 
negatively charged Asca — X observed in this work. 

Moreover, the electronic structure of Asca — Vas has 
been studied by total-energy Green's-function calcula- 
tions treating many-body effects within LSDA-DF^--: and 
shows that the electronic levels allow at most — 1 charging 
for this complex which is in agreement with our results. 

On the experimental side, negatively charged Vas — X 
complexes have been recently detected by positron an- 
nihilation experiments^ between 20—330 K in weakly 
p-type GaAs under arsenic rich condition. The detection 
of Vas — X complexes in these samples was surprising and 
have never been reported before since under these con- 
ditions^S isolated arsenic vacancies arc unlikely to form. 
Vca diffusion to the first neighbor after annealing could 
explain these finding assuming that the metastablc com- 
plex Asca — Vas could be detected at sufficiently low 
temperatures. 

In addition, Asca — Vas structure was proposed by 
Steinegger et al^ to be a potential candidate for EL6 
defect. Measuring the relative concentration of EL6 by 
photon-induced current-transient spectroscopy at room 
temperature, they observed an increase of AsQa concen- 
tration with annealing time by roughly a factor of 2 cor- 
responding to a to decrease of Vca concentration from 
about 10^^ to zero. While the first nearest-neighbor dif- 
fusion cannot play a role in self-diffusion, we suggest that 
it is dominant in these conditions, leading to the trans- 
formation of a large number of Vca into AsQa — Vas com- 
plexes. 



B. Diffusion to second neighbor 

Nevertheless, diffusion to the second neighbor is more 
interesting since it conserves the equilibrium concentra- 
tion of vacancies. The plane-passing paths of Vca simu- 
lated with SIEST-A-RT confirm that the transition state 
does cross the diffusion plane in all charge states. For all 
our simulations the As^^* atom crosses the plane before 
Ga^"'' atom; we did not record any event where As^^^ 
crossed out the plane or where it was located on the plane 
as suggested by Bockstedte and Scheffler— . However, the 
presence of the extra charge affects the diffusion trajec- 
tory by increasing the distance between the moving atom 
(Ga^"'') and the three As^''* dangling bonds. For V^^ 
and Vq^ the charge on As^"*' is so strong that is scatters 
the moving atom and pushes it away from the vacancy 
toward the plane. Consequently, the distance between 
the plane and the atom is lowered as the charge of the 
vacancy increases, it almost vanishes for —3 charge state. 

The diffusion barrier of plane-passing mechanism is 
found to increase moderately as extra electrons are 
added. Recent theoretical works reported a charge- 
dependent migration to the second neighbor by plane- 
passing mechanism for vacancies in SiC and GaN binary 
semiconductors. Bockstedte et ali^ found that the mi- 
gration barrier for diffusion Vc and Vsi in 4C-Si decreases 
when the vacancy charges are increasing progressively. A 
similar trend was observed for Vn in GaNi^ and recently 
for Voa in GaN^. Our most recent results with SIEST- 
A-RT concerning Vas in GaAs^ suggest that the mi- 
gration barrier by plane-passing mechanism decreases by 
increasing progressively vacancy charging. Thus, Vca in 
GaAs shows an opposite trend compared to previously 
studied vacancies. This trend was observed previously 
by Bockstedte and Schefflei>ii as they found a migration 
barrier of 1.7 eV for neutral and 1.9 eV for —3 charge 
state. This behavior cannot be a drawback of SIEST-A- 
RT, since the total path relaxed using CI-NEBM lead to 
the same barrier, it is rather correlated to the electronic 
charge diffusion of Vca in GaAs. 

Experimentally, negatively charged gallium vacancy 
migration on GaAs(llO) surface have been found to be 
stimulated by the tip during STM experiments^ only 
when holes are injected on the surface reducing conse- 
quently the negative charging of the vacancy. This sup- 
ports our results suggesting that the migration barrier is 
lowered as the vacancy becomes less negatively charged. 

C. Diffusion in experimental systems 

Bliss et alMi identified the migration barriers for Vca 
in LT-GaAs by positron annihilation technique. They 
found a migration enthalpy of 1.5±0.3 eV for the vacancy 
diffusion to the second neighbor and 1.1±0.3 eV for dif- 
fusion to the first neighbor. Within the experimental er- 
ror bars, these results agree with our calculated barriers 
summarized in Table UTTI Our results are also consistent 
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with the widely accepted value of 1.7-1.8 eV for GaAs 
vacancy migration in bulk and GaAs-based QW mate- 
j.jg^jg44j4^j4g,4j^^ This value has been recently reported 
for intcr-difFusion in InGaAs/GaAs quantum dotsi^ as 
well. Thus our calculated diffusion barriers are accurate 
enough and could be extended to analyze experimental 
data in bulk and even in nanostructured materials. 

Our study can be useful for interpreting more accu- 
rately experimental data, especially if taking into ac- 
count the fact that most experiments for gallium self- 
diffusioniSiS have been performed in n-type or intrin- 
sic material where negatively charged vacancies are more 
abundant. The most relevant pathway for self-diffusion 
is the plane-passing mechanism since it has the lowest 
barrier for diffusion to the second neighbor, while at suf- 
ficiently low temperatures diffusion to the first neighbor 
could be observed for —1 and —2 charge states. By com- 
bining information coming from our calculated first and 
second neighbor diffusion barrier it becomes possible to 
identify the charge state of the diffusing vacancy in some 
experiments. For example, in the experiment of Bliss et 
alm^, Voa diffuses to the first neighbor by crossing a bar- 
rier of 1.1±0.3 eV, which is close to our calculated barrier 
of 0.84 eV for Vq^ . In addition, measured diffusion bar- 
rier to the second neighbor is 1.5±0.3 eV which is compa- 
rable to 1.7 eV we calculated for —1 charge state. Thus, 
we can state that diffusing vacancies in this experiment 
are more probably singly negative instead of the generally 
accepted triply negative. Recent research works support 
more and more the hypothesis of —1 charge state. This 
was the case, for example, in a recent theoretical study 
comparing simultaneously the relaxations and lifetimes 
obtained from simulation and Doppler spectra of positron 
annihilation experiments, the results confirm the possi- 
bility of —1 charge state for Voa rather that the —3 charge 
stated. 



V. CONCLUSIONS 



The SIEST-A-RT approach has been used to study var- 
ious diffusion pathways for charged gallium vacancies in 
GaAs demonstrating that diffusion in bulk semiconduc- 
tors is a rich and complex phenomenon closely related to 
the charge state. Novel diffusion pathways are identified 
for negatively charged vacancies, showing that diffusion 
process could be non intuitive and more complex than 
initially thought. Voa can diffuse to the second neighbor 
using two different mechanisms in addition to diffusion 
to the first and fourth neighbor. Diffusion to the second 
neighbor is possible by the plane-passing mechanism at 
all charge states as well as by the newly found cluster- 
assisted mechanism which becomes possible only in the 
presence of negatively charged vacancies. In addition, 
gallium vacancy diffusion to the first nearest neighbor is 
possible only for g = 0, — 1, — 2 charges passing by regular 
or distorted Asca ~ AT defects. However this mechanism 
cannot be considered to be dominant for Vca since the 
complete diffusion process was impossible to achieve. Fi- 
nally, the highest diffusion barrier was recorded for the di- 
rect diffusion to the fourth neighbor along the 100 direc- 
tion for negatively charged vacancies. All barriers for the 
migration pathways of Vca are found to increase with in- 
creasing the number of electron, a behaviour opposite to 
what was recently found in the case of vacancy-mediated 
self-diffusion in SiGi2ii^ and GaNii^. 
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